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Galaxies Inside Stromgren Spheres of Luminous Quasars at z > 6: 

Detection of The First Galaxies 

Renyue Cen 1 
ABSTRACT 

The intrinsic Lyct emission lines of normal galaxies before reionization are 
much absorbed by the damping wing of the Gunn-Peterson trough, rendering 
their direct detection nearly impossible, if their intrinsic line widths are less than 
~ lOOkm/s. High redshift luminous quasars prior to the completion of cosmo- 
logical reionization at z ~ 6, on the other hand, are capable of producing large 
HII regions around them (Stromgren spheres) to allow their intrinsic Lya emis- 
sion lines to be transmitted without overwhelming absorption (Cen & Haiman 
2000). We suggest that targeted observations at the Stromgren spheres of known 
luminous quasars at z > 6 would be able to detect Lya emission lines of galaxies 
inside the Stromgren spheres largely unattenuated. A tunable, very narrowband 
filter of ^ ~ 0.1% or a narrowband filter of 4^ ~ 1% with follow-up spectro- 
scopic identifications will be required. Such observations could directly observe 
the sources of cosmological reionization including possibly the Pop III galaxies 
at z = 6 — 20 by JWST. Possible applications include determinations of the 
ionization state of the intergalactic medium, the sizes of the Stromgren spheres, 
the ages of the quasars, the luminosity function of high redshift galaxies and 
its evolution, the spatial distribution of galaxies and its evolution, the biased 
distribution of galaxies around quasars and the anisotropy of quasar emission. 
Observations using Keck-class telescopes may already be made to enable a dif- 
ferentiation between a fully neutral and a 10% neutral intergalactic medium at 
z > 6. 

Subject headings: cosmology: theory — intergalactic medium — reionization 
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1. Introduction 

Recent observations of high redshift quasars from the Sloan Digital Sky Survey (SDSS) 
suggest that the cosmological reionization was completed at z ~ 6 (e.g., Fan et al. 2001; 
Becker et al. 2001, Barkana 2001; Djorgovski, Castro, & Stern 2001; Cen & McDonald 2002; 
Litz et al. 2002), implying a largely opaque universe beyond z ~ 6. In the context of the 
standard cold dark matter model (Bahcall et al. 1999) there is by now a clear consensus that 
galaxies of mass 10 8 — 10 9 M are largely responsible for the final cosmological reionization 
(e.g., Barkana & Loeb 2001; Madau 2002; Wyithe & Loeb 2003a,b; Mackey, Bromm, & 
Hernquist 2003; Cen 2003a,b; Venkatesan, Tumlinson, & Shull 2003; Ciardi, Ferrara, & 
White 2003; Somerville & Livio 2003; Fukugita & Kawasaki 2003; Sokasian et al. 2003). The 
recent Wilkinson Microwave Anisotropy Probe (WMAP) observations (Kogut et al. 2003) 
appear to confirm a much more extended reionization process, with the first reionization 
occurring at a redshift possibly as high as z = 15 — 25 (Cen 2003a). As a result, our 
observational frontier is quickly being pushed towards high redshifts not thought possible 
just a few years ago. It is of great interest to directly detect the sources of reionization and 
the ionization state of the intergalactic medium from z = 6 to z = 30. 

To date the observed Lyman alpha emitters (LAE) at 6 < z < 6.6 appear to be quite 
luminous with star formation rate (SFR) probably in excess of 40M Q yr _1 (Hu et al. 2002; 
Kodaira et al. 2003), which, if assuming a starburst of duration of 5 x 10 7 yrs and star 
formation efficiency of 20%, would give > 1 x 10 10 M Q of stars. The masses of these galaxies 
most likely exceed 1 x 1O 11 M . These observed LAEs are thus much more massive than 
the typical galaxies thought to produce most of the ionizing photons at z > 6 - current 
observations are just scraping the tip of an iceberg. Because of the damping wing of the 
Gunn-Peterson trough cast by the significantly neutral intergalactic medium (IGM) before 
reionization (Miralda-Escude 1998), the intrinsic luminosities (after correcting for possible 
gravitational lensing effects) for these observed LAEs may be still larger. Consistent with 
this, typical LAEs observed at somewhat lower redshift appear to be less luminous by a 
factor of ten or more (Yan, Windhorst, & Cohen 2003; Rhoads et al. 2003; Taniguchi et al. 
2003). For smaller galaxies of mass 10 8 — 10 9 M typically with an intrinsic Lya line width 
narrower than lOOkm/s, the absorption by the damping wing of intervening neutral IGM is 
dramatically more severe, probably rendering them undetectable (Haiman 2002). 

In this Letter we suggest that targeted, very narrowband observations of Lja emission 
at Stromgren spheres of known luminous quasars at z > 6 may provide a unique window to 
directly detect the bulk of the sources for cosmological reionization, because a large fraction 
of galaxies in such regions can transmit their Lyct emission lines largely unattenuated. 
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2. Galaxies Inside Stromgren Spheres of Luminous Quasars 

In the absence of hydrogen recombination with a static density field, the proper radius 
Rs of the Stromgren sphere of a quasar of age tQ is such that the number of hydrogen atoms 
in the sphere equals the total number of ionizing photons produced by the source (Cen & 
Haiman 2000): 
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where (n#) is the mean hydrogen density within Rs, xjgm is the neutral hydrogen fraction 
of the IGM, and N ph is emission rate of ionizing photons from the quasar. For (n#), we have 
adopted the mean IGM density with Q h = 0.047, and N ph = 1.3 x 10 57 s _1 is the observed flux 
of SDSS 1030+0524 (Haiman & Cen 2002). In the realistic case with an evolving cosmological 
density field and a finite hydrogen recombination time, Rs can be computed in detail (Cen 
& Haiman 2000). But Rs estimated in equation (1) gives a reasonable approximation to the 
truth (Haiman & Cen 2002) at redshift z < 8. 

At sufficiently large separations from the Lya line center (^ >2x 10 -8 ) the optical 
depth of the damping wing is (Miralda-Escude 1998) 

r d (AX) - 1.3 x 10 0X)3 ^(^ ( ~6~ ) ( X } ' (2) 

where is baryon density parameter, Y is the helium mass fraction of the primordial gas, 
H and H(z) are the Hubble constant at redshift zero and z, respectively, AA and A are 
the observed wavelength separation from the line center and the observed wavelength of the 
Lya line, respectively. Replacing AA/A by H(z)Ar/c, where Ar is the distance between 
the surface of the Stromgren sphere and the point f inside the Stromgren sphere in question 
along the line of sight and c is the speed of light, and approximating H(z) m H (l + zY^VL 1 ^ , 
we obtain 

r d (r) = l.2x 1GM (—) (— )( — ) , (3) 

where r is the proper quasar-centric distance of the point in question and 9 is the angle 
between the line segment from the quasar to the observer and the line segment from the 
quasar to r. 

The Gunn- Peterson optical depth due to the residual neutral hydrogen inside the Stromgren 
sphere, in the absence of other ionizing sources, can be expressed as (Cen & Haiman 2000) 
r r (r) = 0.34(ff)- 1 /2(_g|_) ( 13x ^f 7 b-1 ) (i±| ) 9 / 2 (^) 2 , due to a very small residual neu- 
tral fraction x r < 10~ 4 . As a result of small x r , the local mini- Stromgren spheres produced 
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by Lya galaxies would render the absorption due to all residual neutral hydrogen inside the 
Stromgren sphere negligible. The total Lya optical depth at r*is hence just r^. From Equa- 
tion (3) we see that, for a luminous quasar like SDSS 1030+0524, Lya lines from galaxies 
located within the central region region of r < IMpc, would suffer a total optical depth of 
r t = 0.22 — 0.36 corresponding 9 from 180° to 0°, i.e., a flux reduction of < 31%, assuming 
x\gm = 1-0. It is evident that there is a significant region within the Stromgren sphere 
where most of the Lya emission galaxies are able to transmit the Lya emission lines largely 
unattenuated. 

There are several potentially very interesting applications from direct observations of 
regions around luminous quasars prior to the completion of cosmological reionization. Be- 
fore turning to the possible applications it is useful to estimate the abundance of possi- 
bly detectable galaxies at z > 6. Assuming that the fraction of all mass having formed 
into halos with masses M h > 10 9 M Q is f g and that starbursts each with a duration of 
t s b has occurred once in each halo in the last half Hubble time, then we find that the 
comoving number density of such galaxies is n g = l-2( 5x *Q7 yrs )(^Q)( 10 j^ I )~ 1 h 3 Mpc~ 3 . 
Assuming a star formation efficiency c* we obtain the star formation rate in each galaxy 
SFR = 0-68(g-|j)( 5x iQTyrs )" 1 ! i 9m q ) M Q /yr. This simple estimate, however, involves two 
currently unknown parameters, c* and t^. 

Perhaps a better way to estimate the abundance and luminosity of LAE is to utilize 
both hierarchical CDM model and LAE survey observations at highest accessible redshift, 
currently at z ~ 5.7. Figure 1 shows the predicted cumulative number density of galaxies at 
z = 5.7 (solid curve) and z = 6.5 (dashed curve) for a cosmological constant dominated cold 
dark matter model with Vt M = 0.29, A = 0.71, H = 70km/s/Mpc, a 8 = 0.80 and n = 1.0. 
The observational data point for LAE shown as the solid dot in Figure 1 is taken from 
Rhoads et al. (2003), who give n(SFR > 5.3 M©^ -1 ) = 5.9 x Kr 5 Mpc -3 at z ~ 5.7. We 
have "normalized" the observational data point to the theoretical prediction at z = 5.7 by 
requiring the same abundance, which then gives an implied mass threshold of the observed 
LAE by Rhoads et al. Clearly, even the LAE observed at z ~ 5.7 by Rhoads et al. are quite 
massive with masses about 1O 12 M , if the adopted cold dark matter model is correct. The 
computed lower bound of star formation rate for the observed LAE at z ~ 5.7 by Rhoads 
et al. of SFR > S.SMQyr^ 1 corresponds to the flux sensitivity of 2 x 10~ 17 erg/cm 2 /sec of 
their LALA survey. If we assume that the Lya luminosity of an LAE is simply proportional 
to its mass, it would require a flux sensitivity of (10~ 18 , 10~ 19 , 10~ 20 )erg/cm 2 /sec to detect 
galaxies of mass of (10 11 , 10 10 , 10 9 ) M with corresponding abundance of (0.01, 0.3, 5.0)Mpc~ 3 
at z — 6.5 (open circles in Figure 1), respectively. 

Between night sky lines there are a few good observational windows at A > 0. 79 /im (cor- 
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Fig. 1. — shows the cumulative number density of galaxies as a function of mass at z — 5.7 
(solid) and z = 6.5 (dashed) in a cosmological constant dominated cold dark matter model, 
consistent with WMAP observations, with the following parameters: Qm — 0.29, A = 0.71, 
H = 70km/s/Mpc, <j 8 = 0.80 and n = 1.0. The solid dot is an observational data point 
from Rhoads et al. (2003). The open circles are three representative galaxy masses with 
approximate star formation rates and Lya line fluxes indicated. 
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responding to z > 5.5 for the Lya line), most notably the one center at z ~ 5.7 and the other 
at z ~ 6.5. However, even within those "good" windows the required integration time to 
detect fainter LAE using Keck telescopes is quite long, t = (flux/10~ 17 erg/cm 2 /sec) -2 hours 
(Gjorgovski 2003, private communications; Fan 2003, private communications), implying 
that approximately a hundred hours is needed on Keck to detect galaxies of mass 1O 11 M , 
even in the absence of significantly attenuation by neutral IGM. The current observational 
limit for Lya flux roughly corresponds to SFR ~ a few M /yrs at z ~ 6 (Hu et al. 2002; Yan, 
Windhorst, & Cohen 2003; Rhoads et al. 2003; Taniguchi et al. 2003). Detecting even lower 
mass LAE may be quite infeasible using Keck telescopes. Furthermore, for ground-based ob- 
servations where the primary background noise in the interested spectral regime comes from 
sky background, it is a very tedious task to search for usable windows between numerous 
strong sky lines to find high redshift galaxies and quasars. Hence, it is likely necessary to 
resort to James Webb Space Telescope (JWST) for detecting the bulk of galaxies (of mass 
M h ~ 10 9 M Q ) at z > 6. As pointed out by Oh (1999), JWST should be able to detect the 
H a line (and possibly H^ and other Balmer lines) of galaxies of mass Mh ~ 10 9 M Q at z > 6. 
A joint detection in Lya and H a would not only firmly determine the redshift of each galaxy 
but also shed very useful light on important physical processes of each galaxy including dust 
content. In addition, the ratio of Lya to H a lines would provide a cross-check of the damping 
wing modulated spatial distribution of LAE or a tool to break possible degeneracies in the 
spatial distribution of LAE from Lya line alone. 

We note that three factors could boost the effectiveness of star formation at high redshift 
in terms of Lya emission. First, the gas at high redshift may be much metal poor; for 
example, taking the same Salpeter IMF would give about a factor of 2 more ionizing photons 
hence Lya photons for stars with Z m 10~ 2 Z Q compared to stars with solar metallicity. 
Second, if the metal poor gas at high redshift forms relatively more massive stars, the 
efficiency of producing ionizing photons would be higher. Third, the dust content in the 
high-z metal-poor galaxies should be lower, possibly allowing for relatively more Lya photons 
to escape. These factors may render high redshift LAE more luminous that their low redshift 
counterparts at a star formation rate. 



3. Applications 

3.1. Probing the Ionization State of the IGM and the Size of the Stromgren 

Sphere 

Figure (2) shows the ratio of the number of galaxies at a distance AA/A (with width 
S^) from the quasar on the blue side over the number of galaxies at the same distance from 
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Fig. 2. — shows the ratio of the number of galaxies at a distance AA/A (with width 5^-) 
from the quasar on the blue side over the number of galaxies at the same distance from 
the quasar on the red side within a transverse radius of r p = 2Mpc proper, as a function 
of AA/A. Three panels correspond to three flux sensitivities of lCT 18 erg/cm 2 /sec (Panel 
a), 10~ 19 erg/cm 2 /sec (Panel b) and lCT 20 erg/cm 2 /sec (Panel c), which are then translated 
into the threshold mass of detectable galaxies using Figure I. In each panel two sets of 
cases are shown. Case (1): R s = 5Mpc proper with xiqm = (1.0,0.1,0.01), respectively, 
shown as solid (green), dotted (red) and dashed (blue). Case (2): R s = 3Mpc proper with 
Xigm = (1.0,0.1,0.01). Each set asymptotically curves downward towards the radius of the 
respective Stromgren sphere indicated by the vertical tick sitting on the bottom x-axis. The 
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the quasar on the red side within a quasar-centric transverse radius of r p = 2Mpc proper, as 
a function of AA/A. Three panels correspond to three flux sensitivities of 10 _18 erg/cm 2 /sec 
(Panel a), 10~ 19 erg/cm 2 /sec (Panel b) and 10 _20 erg/cm 2 /sec (Panel c), which are then trans- 
lated into the threshold masses of detectable galaxies using Figure 1 to infer the correspond- 
ing abundances. In each panel two sets of cases are shown. Case (1): R s — 5Mpc proper 
with xigm = (1-0, 0.1, 0.01), respectively, shown as thick solid, thick dotted and thick dashed 
curves. Case (2): R s = 3Mpc proper with Xiqm = (1.0,0.1,0.01), respectively, shown as 
thin solid, thin dotted and thin dashed curves. We choose to plot the ratio instead of the 
direct number density to avoid ambiguity due to a possibly non-uniform, quasi spherical 
distribution of galaxies around quasars. Galaxies will still likely be clumpy and clustered on 
small scales but a large volume considered would smooth out most of the fluctuations. 

Let us now estimate the accuracy of the data points. The comoving volume enclosed 
by a cylindrical region within a transverse radius r p = 2 Mpc and a line-of-sight distance 
of AA/A = 0.01 (from AA/A = -0.005 to AA/A = 0.005) is 2.1 x 10 4 Mpc 3 comoving at 
z = 6.28, which is then multiplied by the number density of galaxies using the respective 
n from Figure 1 to yield the corresponding total number of galaxies above the respective 
mass threshold. Five data points are used for the indicated AA/A range. The errorbars are 
shown in each panel for the case with R s = 3Mpc proper. From Panel (a) of Figure 2 it 
is clear that xiqm = 1 and xiqm = 0.1 are distinguishable at 2 — 3a level even for galaxies 
with mass M > 10 11 M (Panel a), detectable by Keck telescope with an integration time 
of order of 100 hours for each galaxy (assuming one spectrum per observation) with total 
~ 210 galaxies. This might be a costly project but the potential payoff is attractive. Cases 
with xigm = 0.1 and x\gm = 0.01 can only be marginally differentiated at this sensitivity 
level. At the same time, it is evident that the models with different R s can be more easily 
distinguished, allowing for a precise measurement of R s . In addition, a precise measurement 
of R s , in conjunction with the spatial information (see §3.2), may allow a determination of 
the age of the quasar. With JWST, if one can reach flux levels 10~ 20 — 10~ 19 erg/cm 2 /sec, 
very precise measurements of xigm can be made, as can be seen from Panels (b,c) in Figure 
2. 

We note that AA/A = 0.005 corresponds to a line-of-sight comoving distance from the 
quasar of 16Mpc at z — 6.28. The gravitational influence of the quasar host galaxy at 
that distance is likely to be negligible, alleviating some possible complications including 
infall. In addition, the typical peculiar velocities of galaxies at high redshift is likely be 
< 200km/s, which would correspond to a possible displacement of AA/A = 6 x 10~ 4 , too 
small to significantly matter. 

Finally, we note that, in the case of a lone LAE (not sitting inside a large Stromgren 
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sphere of a luminous quasar), an LAE estimated to have an SFR about an order of magnitude 
higher than the solid dot in Figure 1 (Hu et al. 2002) can only transmit a small fraction of 
Lja flux due to attenuation by the neutral IGM (Haiman 2002). Lower mass LAE are further 
devastated by their lower intrinsic line widths in addition to smaller Stromgren spheres. As 
Haiman (2002) pointed out, no Lya line of an LAE narrower than lOOkm/s can be observed. 
Therefore, it is estimated that no LAE galaxies less massive than 10 10 — 10 11 M will be 
detectable in a random field outside a large Stromgren sphere of a luminous source. 

3.2. Galaxy Luminosity Function and Spatial Distribution 

Fitting the data to obtain R s and xiqm (§3-1) would yield precise information on r^, 
which in turn would allow a determination of a region on the red side of the quasar where 
attenuation due to intervening neutral IGM is negligible. Consequently, one can construct 
an unbiased luminosity function given the likely availability of up to thousands of galaxies in 
that region, which would further allow for corrections to be made for "missing" galaxies in 
other regions where is appreciable. Moreover, the spatial distribution of galaxies in this 
region and others in other regions would allow an accurate determination of the clustering 
properties of galaxies up to a scale of 10 — 30Mpc comoving. 

Given the determined luminosity function and £igm, along with follow-up spectroscopic 
observations of the galaxies detected which would produce the line widths for these galax- 
ies, we could self-consistently determine a fraction of these galaxies that would be directly 
observed even sitting in an average region and their possible observable features. 

3.3. Biased Distribution Around Quasars 

Taking the farthest red region as an unbiased region, this would then allow a determina- 
tion of luminosity-dependent biased distributions of galaxies, down to the lowest luminosity 
of galaxies in the farthest red region. This may shed useful light on the formation of quasars 
and their host galaxies. Whether or not the farthest red region is an unbiased, "average" 
region can be checked using the same data, where one may expect to see a diminished radial 
gradient. 
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3.4. Probing Anisotropy of Quasar Emission 

The three-dimensional information of the distribution of galaxies around a quasar makes 
it possible to probe the angular distribution of emission from the quasar. It should be rela- 
tively straightforward to build self-consistent models to account for possible outcomes from 
observations. For example, if galaxy abundance drops off rapidly at a transverse distance 
significantly smaller than the inferred R s , it would indicate that the quasar emission is bi- 
polar, whose angles can be inferred. If, in addition, there is a strong suppression of galaxies 
on the red side of the quasar even at small r p , it would imply that the quasar emission is 
largely beamed towards us. We will not attempt to make more detailed models here. If 
the standard unified model for AGN/quasars (e.g., Antonucci 1993) is correct, the ensuing 
consequence on the distribution of LAEs around quasars will clearly be detectable. This 
method may provide the first and accurate way to test this picture. 

In summary, we show that, with a tunable, very narrowband filter of AA/A ~ 0.1% or a 
narrowband filter of ^ ~ 1% to survey a luminous quasar Stromgren sphere with follow-up 
spectroscopic identifications of the Lya emission lines of individual galaxies inside the quasar 
Stromgren sphere, we could obtain a wealth of information for the high redshift galaxies that 
are responsible for reionizing the universe at z ~ 6. This same technique may be applicable 
to quasars at any high redshift. At higher redshift, galaxies are somewhat less massive, but 
they may have a more top-heavy IMF, possibly a higher star formation efficiency, a shorter 
starburst duration and much lower dust content. These positive, compensatory effects could 
possibly render high redshift galaxies with mass 10 7 — 10 8 M intrinsically as luminous as 
galaxies of mass 10 8 — 10 9 M at z — 6. We expect that JWST should be able to detect 
these galaxies at z > 20 (Eisenstein 2003), provided that a tunable, narrowband filter at 
appropriate wavelengths is available. 

4. Conclusions 

We show that targeted observations at the Stromgren spheres of known luminous quasars 
with tunable very narrowband filters of width ~ 0.1% or a narrowband filter of ^ ~ 1% 
to survey a luminous quasar Stromgren sphere with follow-up spectroscopic identifications 
of the Lja emission lines of individual galaxies would provide a direct way to observe the 
Lya emission lines of the bulk of galaxies at z > 6 in the mass range 10 8 — 10 9 M Q that are 
believed to be responsible for the second cosmological reionization. At present it appears 
that ground-based observations, such as by Keck, can only detect a relatively small fraction, 
high end tail of galaxies with mass > 1O 11 M . But even such observations may already 
be able to differentiate between a fully neutral and a 10% neutral intergalactic medium at 
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z > 6. JWST should provide the unique opportunity to detect the bulk of galaxies at high 
redshift. The likely gain in sensitivity with a tunable, very narrowband filter may allow 
JWST to detect the Pop III galaxies of mass 10 7 - 10 8 M Q with a top-heavy IMF at z > 20. 
In this case, the key is to search for luminous enough quasars at these very high redshifts, 
perhaps requiring an infrared survey comparable to SDSS. This is presently technologically 
challenging but the potential scientific payoff is immense. 

We show that such observations will allow for accurate determinations of the ionization 
state of the IGM, the sizes of the Stromgren spheres, the ages of the quasars, the luminosity 
function of high redshift galaxies and its evolution, the spatial distribution of galaxies and 
its evolution, the biased distribution of galaxies around quasars and the anisotropy of quasar 
emission. 

If the universe was indeed reionized twice (Cen 2003a), then there will be a window 
of opportunity somewhere at z = 10 — 25, where galaxies can transmit Lja lines (as well 
as optical-UV continuum and other emission lines) unabsorbed; they should be directly 
observable by JWST even in the absence of luminous quasars. A signature continuum cutoff 
of the spectra of these galaxies shortward of 1216(1 + z rec ^) 'k would pin down the redshift 
of the second cosmological recombination epoch z reCt 2 and then a rapid thinning out of the 
number density of Lya galaxies at some redshift z > z reCt2 would be indicative of the first 
reionization epoch. 

This research is supported in part by NSF grant AST-0206299. I thank George Djor- 
govski, Xiaohui Fan and Rogier Windhorst for useful comments and discussions, Michael 
Strauss for pointers to useful astronomy websites, and the referee for useful suggestions 
which significantly strengthen the paper. 
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